Large amplitude behavior of the bulk viscosity of dense matter 
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We study the bulk viscosity of dense matter, taking into account non-linear effects which arise 
in the large amplitude "supra-thermal" region where the deviation /xa of the chemical potentials 
from chemical equilibrium fulfills /ia > T. This regime is relevant to unstable modes such as r- 
modes, which grow in amplitude until saturated by non-linear effects. We study the damping due to 
direct and modified Urea processes in hadronic matter, and due to nonleptonic weak interactions in 
strange quark matter. We give general results valid for an arbitrary equation of state of dense matter 
and find that the viscosity can be strongly enhanced by supra-thermal effects. Our study confirms 
previous results on quark matter and shows that the non-linear enhancement is even stronger in the 
case of hadronic matter. Our results can be applied to calculations of the r-mode-induced spin-down 
of fast-rotating neutron stars, where the spin-down time will depend on the saturation amplitude of 
the r-mode. 



I. INTRODUCTION 

Compact stars are the only known objects that contain 
equilibrated matter that is compressed beyond nuclear 
density, making them a valuable laboratory for the study 
of the structure of matter under extreme conditions. In 
addition to hadronic matter they may also contain new 
forms of matter that involve deconfined quarks [IrH] • In 
contrast to the static properties of compact stars which 
only depend on the equation of state of matter [5J, dy- 
namic properties also depend on the low energy degrees of 
freedom and thereby might be able to discriminate more 
efficiently between different forms of strongly interacting 
matter. One of the dynamic properties of dense matter 
is viscosity, which determines the damping of mechanical 
perturbations, and a particularly important application 
is to the damping of r-mode oscillations of compact stars 
[tjlllOj. which, at sufficiently low viscosity and high ro- 
tation rate, are unstable and can cause rapid spin-down 
of the star |11| via gravitational radiation. Since the r- 
mode is unstable, its exponential growth must eventually 
be stopped by some non-linear mechanism. Finding the 
relevant mechanism is important because it determines 
the amplitude at which the r-mode saturates, and hence 
the rate at which it spins down the star. Previously sug- 
gested mechanisms include mode coupling and the trans- 
formation of the r-mode energy into differential rotation 
P2HH], and friction between different layers of the star 
such as "surface rubbing" at the crust [TB]. Because of 
the complexity of the problem, these mechanisms have 
to rely on approximations that are not always well con- 
trolled. 

In this paper we consider an alternative mechanism 
which does not involve additional physics, but is al- 
ready present in a quasi-static hydrodynamic descrip- 
tion. At low amplitudes the bulk viscosity is amplitude- 
independent, but since the r-mode is unstable its am- 
plitude grows, and unless stopped by other mechanisms 
will quickly enter the "supra-thermal" regime where the 
bulk viscosity grows with amplitude, and may become 
large enough to stop the growth of the mode [TTl H8] , 
The supra-thermal regime is characterized by /ia ^ T, 



where T is the temperature and /iA is the amplitude of 
the oscillations in the chemical potential of the quantity 
whose re-equilibration causes the viscous damping. We 
will study the microscopic part of the problem via a com- 
prehensive analysis of the bulk viscosity of dense matter. 
We leave the astrophysical aspects for future work. Since 
the precise phase structure and the equation of state of 
matter at high density is still unknown, we keep the de- 
pendence on those parameters as explicit as possible, and 
as well as numeric results we provide analytic approxima- 
tions which prove to be surprisingly accurate. This allows 
us to obtain general results for the bulk viscosity valid 
for many different phases of matter, and enables us to 
estimate the involved uncertainties. We study in detail 
the cases of equilibrated npe-matter and strange quark 
matter, and consider both modified and direct Urea pro- 
cesses in the hadronic case. Yet, our general expressions 
can be applied to other equations of state and entirely 
different forms of strongly interacting matter. 



II. BULK VISCOSITY OF DENSE MATTER 

The bulk viscosity of a given form of matter is a mea- 
sure of the energy dissipated when it is subjected to an 
oscillating cycle of compression and rarefaction. Bulk vis- 
cosity is known to be the dominant source for the damp- 
ing of r-mode oscillations at high temperatures and low 
amplitudes. Consequently, bulk viscosity has been com- 
puted for many forms of dense matter [171 [T^H57] . Nearly 
all these studies restricted themselves to the sub-thermal 
regime /ia <C 2ttT. However, as noted above, the astro- 
physically interesting scenario is one where the r-mode 
is unstable, and so the supra-thermal bulk viscosity may 
well become relevant. The influence of the supra-thermal 
regime has been studied numerically in |17| for the case of 
strange quark matter. This analysis showed that for large 
amplitude oscillations the viscosity can increase by orders 
of magnitude. Yet, because of their qualitatively differ- 
ent low energy degrees of freedom and weak-interaction 
equilibration channels, other forms of matter could show 
different behavior. 



In the following we will derive the non-linear equa- 
tions that determine the bulk viscosity due to weak in- 
teractions that interconvert the fermionic species that are 
present 1 . We will then solve it for arbitrary amplitudes. 
We focus on weak interactions because their equilibration 
rate is comparable to typical compact star oscillation fre- 
quencies: strong interactions make a negligible contribu- 
tion at these frequencies because their equilibration rate 
is much too fast. 

The bulk viscosity of a given form of matter is defined 
by the response of the system to an oscillating compres- 
sion and rarefaction. This corresponds to an oscillation 
in the densities of all exactly conserved quantities. We 
will assume that there is at least one such quantity whose 
density we call n*. In compact stars it is typically the 
baryon number density. We will study the energy dis- 
sipated as a result of a small harmonic oscillation Sn* 
around its equilibrium value n* 

n* (r, t) = n*(r) + Sn* (f, t) = n*(r) + An* (f) sin f ) , 

(1) 
where the amplitude of the oscillation is An*, and we 
assume An* -C h*. The dissipated energy per volume 
due to the oscillation is given by 



de 
dt 



<(v-«)' 



(2) 



where v is the local velocity of the fluid of the conserved 
quantity and the continuity equation for its particle num- 
ber n* reads 



we can express the dissipated energy per volume averaged 
over one time period in terms of the induced pressure 
oscillation 



1 f T p dn* 
n* dt 



dt. 



(7) 



To calculate the bulk viscosity we must calculate p(t) . 
We will assume that the bulk viscosity arises from beta- 
equilibration of fermionic species. We further assume 
that, in the absence of weak interactions, there would be 
s conserved species, and that there is a single channel 
of weak interactions that can perform interconversion of 
species, leaving s — 1 exactly conserved fermion-numbcr 
charges 2 . We defer discussion of the general situation 
of several coupled channels to future studies. Subtract- 
ing the chemical potentials of the final state particles in 
the relevant weak channel from those of the initial state 
particles, we obtain the difference 



/'-A 



= E^-E 



/';/■ 



(8) 



/ 



which is the quantity that is driven out of equilib- 
rium by the driving density fluctuation, and whose 
re-equilibration leads to bulk viscosity. The quasi- 
equilibrium state can generally be described in terms of 
the driving density n* and the ratio x = n\/n* where n\ 
is the density of one of the particle species whose number 
is changed by the equilibration process. For small oscil- 
lation amplitudes An,/fi t <C 1 the pressure can then be 
expanded around its equilibrium value p — p (n*) 



dn 

~di 



- + V-(n*u)=0. 



(3) 



In the case that density varies slowly enough so that den- 
sity gradients can be neglected, and using An*/n* -C 1, 
averaging over a whole oscillation period r = 2ir/cu gives 
the bulk viscosity as 



(An*r 



(4) 



Using the relationship between fluctuations in volume 
and fluctuations of a conserved quantity, 



p = p + 



dp 

dn* 



dp 
Ox 



Sx , 



(9) 



where Sx is the deviation of x from its beta-equilibrium 
value. The t- independent part p as well as the term pro- 
portional to the driving density fluctuation Sn* do not 
contribute to the viscosity integral. The remaining sus- 
ceptibility can be rewritten 



dp 
dx 



d"A 
dn* 



(10) 



Because of weak interactions, x depends on time, 



dn* 



dV_ 
' V 



and the mechanical work done by a volume change 

P 



de 



V 



-dV 



(5) 



(6) 



Sx(t) = I —dt' . 
K ' Jo dt' 



(11) 
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1 nl [ T dfi. 



n An* J dn* J dt 



dx 



- f dt' cos (uit) dt . (12) 



1 We do not study bulk viscosity arising from the interconversion of 
bosons, such as the light mesons that occur in color-flavor-locked 
phases. 



2 In this counting we exclude fermions like neutrinos, which escape 
from compact stars and so are effectively not conserved. 



We want to point out already at this point that, in con- 
trast to the harmonic driving density oscillation 5n* with 
amplitude An,, the induced chemical potential fluctua- 
tion Sua around the vanishing equilibrium value can have 
a more complicated anharmonic form. 
The fluctuations of the density ratio can be obtained from 
an analogous expansion of the chemical potential fluctu- 
ation 



Sfl 



A 



dn* 



9/j.a 

dx 



Sx 



which yields a linear equation relating ua and Sx 

dx 



du A n An* 

— : — = Cw— — cos (ut) 
at n* 



Bn» 



dt ' 



with the susceptibilities 
dfiA 



C = 



dn* 



B = 



1 dn L 



dx 



(13) 



(14) 



(15) 



Using (|14| we obtain the bulk viscosity in terms of the 

-> 

HA(t) cos(u>t)dt . (16) 



chemical potential fluctuation 

1 n* c r 

/(I 



c = — 



7T ArU B 



In terms of a Fourier expansion of the periodic chemical 
potential fluctuation 



Ma(*) 



oo 



(a n sin (nut) + b n cos(nut)) (17) 



we see that the only component of ^a ( l ) that contributes 
to the viscosity is the component of the fundamental 
Fourier mode that lags the driving volume oscillation by 
a phase of tt/2. This suggests that a truncated Fourier 
ansatz may provide a reliable approximation for the vis- 
cosity: we will explore this idea in Sec. |III| 

To obtain the temperature and amplitude dependence 
of the bulk viscosity, we now discuss the general form of 
the beta equilibration rate. We define the net equilibra- 
tion rate 



roo = r ( ^ -i* 1 -) 



dx 



(18) 



where we use the convention that T^> is the rate for 
the process where n\ is decreased, and r^) is the rate 
for the inverse process. We study equilibration processes 
where the net rate takes the general form 



N / 2 



(19) 



where N is the highest power of ua arising in the rate. 
In terms of dimensionless variables 



(p = ut 



M<f) 



MA (t) 
T 



(20) 



the differential equation for the chemical fluctuation 



eq. ( 14 ) can be written as 



dA 
dip 



= d cos (ip) 




(21) 



with the prefactors of the driving and feedback term 
given by 



d = 



C An 
T n* 



/ = 



BTT K 



(22) 



Note that the feedback term involves both linear and 
non-linear parts which are controlled by a single param- 
eter / and that its particular form is determined by the 
constants \j which parametrize the particular weak rate. 
The viscosity is then finally given by 



c = 



TC 



iB Ar 



271 



-4 (ip; d, f) cos (ip) dip , 



(23) 



where A is the periodic solution to eq. (21 1 



Before we discuss the general solution of these equations 
in detail, let us consider its asymptotic limits. 
Sub-thermal limit 

In the limit /xa <C T corresponding to A ^C 1 the non- 
linear terms can be neglected 



d 

dip 



f ) A = dcos(ip) . 



(24) 



Since this equation is linear, the fluctuation A must be 
harmonic and only the n = 1 term in the Fourier ansatz 
eq. (17 1 is present. Inserting this ansatz yields the solu- 



tion for the required Fourier coefficient 



6i 



df 



1 + /2 



(25) 



Inserted in eq. ( 23 1 this yields the general sub-thermal 



result, denoted by a superscript <, for the bulk viscosity 
of an arbitrary form of matter which shows the charac- 
teristic resonant form 



C 2 TT H 



(BYT K ) 2 



imax n 



2uBTT K 



(BTT 



re\2 



(26) 



As long as the combination of susceptibilities C 2 jB does 
not vary too quickly with temperature, the sub-thermal 
viscosity has a maximum 



■>max 



c 2 

2uB 



at J. ry 



u 
TB 



(27) 



Supra-Thermal limit 

The opposite, suprathermal limit, /ia 3> T, corresponds 
to A ^> 1. Since the feedback term in the differential 
equation is restraining, this limit can only be reached in 



the limit of large driving terms d ^> 1. In this case only 



the largest power of A is relevant and eq. (21 I reduces to 

i 

An* 



= dcos((p)-XNfA 



2N+1 



A 



(28) 



The viscosity scales correspondingly in this limit as 
'An, 



C 



2N 
' 2N+1 



(29) 



and decreases at very large amplitudes. 
General solution 

After these limiting cases we will discuss the qualitative 
aspects of the general solution eq. ( 23 1 . Due to the non- 
linearity of the differential equation (21 1 this requires a 
numeric solution. Yet, for each weak channel, character- 
ized by the constants Xji such a solution as a function 
of the two independent variables d and / has to be per- 
formed only once and is then valid for any equation of 
state and includes the complete dependence on the un- 
derlying parameters in eq. ( 22 1 . The qualitative behav- 



ior of the solution as a function of the two independent 
parameters d and / is shown for hadronic matter with 
modified Urea process in fig. [T] Turning up the feedback 
term at fixed driving term increases the phase shift of the 
waveform from to 7r/2 and at the same time decreases 
the amplitude, but the waveform stays harmonic. In con- 
trast, turning up the driving term at fixed feedback in- 
creases the amplitude towards the supra-thermal regime 
A > 1 and the waveform becomes increasingly anhar- 
monic. Recall, however, that only the phase shifted har- 
monic component in the Fourier expansion contributes 
to the viscosity eq. ( 23 1 . 



Motivated by the above expression eq. (27 1 for the 



maximum in the sub-thermal regime the general result 
can be written in the form 

( = C< ax 2(dJ) = ^2(dJ) (30) 

where the dimensionless function X that includes the non- 
trivial parameter dependence is given by 



l(d,f) = 



ltd 



2tt 



A(tp; d, f) cos((p)d<p (31) 



The expression I can then be tabulated as a function of 
the independent parameters d and /. We believe that 
presenting our results in this form will make them easier 
to apply to calculations of r-mode damping, where the 
complete parameter dependence is required. The com- 
putation of the damping time of the mode involves an 
integral over the star of an expression that involves the 
bulk viscosity (e.g., [TT]) which varies throughout the 
star because of its dependence on the amplitude of the 
mode and the susceptibilities, both of which are position- 
dependent. The function X{d, /) encapsulates the de- 
pendence of the bulk viscosity on the position-dependent 



parameters, allowing straightforward evaluation of the 
damping time integral. 

The function X(d, /) is shown in fig. El for two exam- 
ples: a model of strange quark matter and a model of 
hadronic matter; details of the models are discussed be- 
low. We see that the function has the same qualitative 
form in both cases. It has a global maximum value of 
1, reached in the sub-thermal limit and a line of local 
maxima along a parabola in the d- /-plane. Thus the 
maximum value ( 27 1 of the sub-thermal viscosity is also 



the maximum in the general case and depends only on 
the equation of state, the density and the frequency but 
is independent of the weak rate. The weak rate influ- 
ences, however, at what temperatures and amplitudes 



the local maxima are reached. As seen from eq. (22), the 



parameter d is directly proportional to the amplitude, so 
that at moderate feedback an amplitude increase does ini- 
tially not affect the viscosity at all, corresponding to the 
amplitude-independent sub-thermal result. But once the 
amplitude becomes sufficiently large we enter the supra- 
thermal regime and the viscosity increases strongly by 
orders of magnitude until it reaches its maximum. The 
size of the amplitude A is denoted in fig. [2] by the dark- 
ness of shading of the surface. This qualitative behavior 
has already been observed in [T7] but we find that at even 
higher amplitudes the viscosity decreases again accord- 
ing to the limiting behavior eq. ( p9| . In contrast, at large 
feedback the viscosity becomes basically amplitude inde- 
pendent over the relevant parameter range as described 
by the sub-thermal result. 

Let us now discuss the dependence of the viscosity on 
the underlying parameters in eq. (22). An amplitude 



increase (keeping all other variables fixed) results in a 
linear increase in the variable d as shown by the dashed 
(blue online) curves in fig. k] An increase in temperature 
changes the viscosity along a line shown by the solid (red 
online) curves. In order to assess the frequency and am- 
plitude dependence we must take into account the pref- 
actor in eq. ( |30| . This prefactor, given by the maximum 
viscosity in the subthermal regime, is shown in fig. [3] for 
the hadronic model of fig. [2Ta). It exhibits a monotonic 
increase with density and inverse angular frequency. An 
increase in angular frequency therefore changes the vis- 
cosity via a change of X towards the negative /-direction 
and furthermore via the overall prefactor featuring an ad- 
ditional 1/w dependence. A density increase has an even 
more indirect impact since it depends on the detailed 
form of the susceptibilities C (n*) and B (n,*) which like- 
wise arise in the prefactors of the viscosity. These depen- 
dencies will be studied in more detail below. 



III. STRANGE QUARK MATTER 

A. General features 

It has been suggested that the true ground state of 
matter at high densities may be strange quark matter 
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Figure 1: Waveform A(4>) = Ma (wt) /T for different values of the two independent parameters. We show only the positive 
half-wave, on a logarithmic scale. Left panel: Fixed driving term d — 1, with varying feedback term / = 0.001, 0.01, • ■ ■ , 1000. 
As / rises, the phase lag increases from zero towards 7r/2, but at the same time the amplitude decreases as 1//. Right panel: 
Fixed feedback term / = 1, with varying driving term d — 0.001, 0.01, • ■ • , 1000. As d rises, the phase lag rises from 7r/4 to n/2 
and the waveform becomes increasingly anharmonic, approaching a square wave in the limit. 



(a) Hadronic matter, modified Urea process 



(b) Quark matter, non-leptonic process 




Figure 2: The function X arising in the general solution eq. (231 for two models of dense matter. Left panel: hadronic matter 



with modified Urea equilibration. Right panel: quark matter with the non-leptonic equilibration process eq. (321. The function 
has a global maximum of 1 reached asymptotically for d — ► 0, / = 1 and a line of slowly decreasing local maxima along a 
parabola in the d-f plane. The shading of the surface denotes the size of the amplitude A so that dark shades of grey represent 
the supra-thermal regime. Eq. ( |22[ ) relates d and / to underlying physical parameters such as temperature T and amplitude. 
An amplitude increase (keeping all other variables fixed) results in a linear increase in the variable d as shown by the dashed 
(blue online) curves. An increase in temperature changes the viscosity along a line shown by the solid (red online) curves. 
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Fig ure 3: The maximum viscosity in the supra-thermal limit 
eq. ( 43 1 of hadronic matter (upper surface) and a hadronic gas 
(lower surface) as a function of baryon density and angular 
frequency. This represents also the analytic prefactor of the 
general expression in eq. «23p. The corresponding plot for 
strange quark matter eq. ( |49[ ) would be trivial since it does 
not depend on the density to leading order and only shows 
the analytic 1/w dependence. 



O [3] consisting of u, d and s quarks. In that case self- 
bound strange stars could exist. In this section we ap- 
ply the results obtained above to strange quark matter, 
which is both an interesting physical scenario and a useful 
introductory example in which one can make illuminat- 
ing simplifications which are not possible for the case of 
hadronic matter discussed in the next section. The bulk 
viscosities of various forms of strange quark matter have 
previously been analyzed |17l |2"9^36 , and the influence 
of strong magnetic fields has recently been discussed [35] . 
The dominant channel in unpaired strange quark mat- 
ter is the non-leptonic flavor changing process 



d + u -f-> s 4- u 



(32) 



The corresponding quark Urea processes, which involve 
leptons, are parametrically suppressed in the ratio T//i q . 
The conserved quantity that tracks the driving oscilla- 
tion can be chosen as the baryon number, with density 
n = i {n s + rid + n u ). The equilibrating chemical poten- 
tial eq. ffity carries in this case the quantum numbers of 
neutral if-mesons and is therefore denoted by 



Hk = ii s 



fJ-d 



(33) 



The rate of the non-leptonic process eq. ( 32 ) is given by 



ro-o = -i 6 - 

q 5tt 5 



G% sin 2 6» c cos 2 9 c ri d n K (4ir 2 T 2 



■&) 



(34) 

From this expression one can directly obtain the equili- 
bration rate parameters T, k and Xi m the parameteri- 
zation eq. (19 1. Their values are given in the first row of 
table |U 



B. Analytic approximation 

Since in this case only cubic non-linearities arise it is 
possible to obtain an approximate analytic solution to the 
non-linear equation (pTj) . Taking into account the above 
observation that only the leading Fourier coefficient in 
the expansion of the chemical potential oscillation con- 
tributes to the bulk viscosity it is natural to seek such a 
solution via a Fourier ansatz up to a given order O 



A(t) 



o 

E 

H,=-0 



J\- n 



(35) 



where the complex form is used to simplify the compu- 
tation. In principle, the amplitude of the leading Fourier 
mode will depend on the truncation order O, but analyt- 
ically solving eq. (21 1 via a computer algebra system to 



order O — 2we find that the coefficients A±2 vanish iden- 
tically. Correspondingly anharmonicities do not directly 
contribute to the viscosity and are even absent to next to 
leading order so that we can restrict our analysis to the 
leading order = 1. Although such a parameterization 
neglects any anharmonicities it properly captures both 
the amplitude and the phase shift of the oscillation even 
in the large amplitude regime. Due to the reality of the 
solution there is only one independent complex Fourier 
exponent determined by a non-linear algebraic equation. 
In the case of quark matter where Xi — f° r * > 1 and 
only the leading non- linear term \ (fi^/T) is present an 
analytic solution of this equation is possible. In this case 
we can decompose the amplitude into real and imaginary 
parts A\ — A R + iAi, obeying coupled equations 



/A i? ,(l + 3x(^, + A?))+A / = 


= -d/2 


(36) 


fA I (l + 3 x (A 2 R + AJ))-A R -- 


= 


(37) 



Note that an analytic solution is only possible because 
the quark matter equations are cubic; other forms of mat- 



ter with higher order non-linearities in eq. (19 1 require a 
numeric solution. The above system of algebraic equa- 
tions has a lengthy analytic solution which we refrain 
from giving here because, as we will see below, it can be 
very accurately approximated by a much simpler expres- 
sion ( 45 ) constructed from a combination of the solutions 



in the sub-thermal and supra-thermal regimes. Therefore 
we now concentrate on the supra-thermal case, denoted 
by the index >, where the temperature-dependent term 
can be neglected, 

A > (if) = 2A R cos ((f) - 2Ai sin (ip) 

M( (q{z?-lf , g(z) 2 1 . \ 

= -y [ vs zq ( Z ) 2 cos{ip) + "^r smM J ' 

(38) 
where the dimensionless quantity z is defined by 



9V3 



Xd 2 f = 



9V3 X fC 2 B K _ 2 /Arc 



T^Ai" 



UJ 



(39) 
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Table I: Weak interaction parameters describing the considered damping process. Here \i q is the quark chemical potential, n is 
the baryon density, no nuclear saturation density and x the proton fraction. 



and 



i 

q(z) = (Vz 2 +i- z y 



(40) 



Eq. ( 23 I then yields the approximate analytic result for 



the bulk viscosity in the supra-thermal regime 



, 2 C 2 J 9V3 X T C 2 B TK _ 2 



C W 3^B«\ 



An, 



(41) 



in terms of the dimensionless function 
9 



Hz) 



-Lz 



sjz 2 + l-z) + (y / z 2 + l + 



-2 



(42) 
This function has a maximum at z max = 3v3. Since 
h(z max ) = 3v3/4; the corresponding maximum value of 
the viscosity is 



2 C 2 ^ C 2 

3^3 Buj {Zmax) ~ 2Buj 



(max = 7r^-rT-. h ( z max) = ^7^ (43) 

which strikingly is the same expression as in the sub- 



thermal limit eq. ( 27 1 . Correspondingly the bulk viscos- 



ity has a universal upper bound C, max that is independent 
of the particular weak damping process. It is directly 
proportional to the oscillation period with a coefficient 
that only depends on the response of the strongly inter- 
acting matter. However, the corresponding temperature 
( 27 1 and amplitude 



An* 
n* 



/ 8w 

3xTT K ~ 2 C 2 B 



(44) 



at which this maximum is reached both depend on the 
weak rate. 

Knowing the upper bound (max and the functional be- 
havior in the extreme sub-thermal and supra-thermal 
limits, allows us to give a simple parameterization of 
the full function for all temperatures and amplitudes. 
We construct a weighted sum of the analytic results in 
the sub-thermal eq. ( 26 I and the supra-thermal regime 
eq. m, 




,wi 



YBj J 3V3u 2 + Y 2 B 2 T 



(45) 



9V3x YBC 2 T K 



8 
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Studies of the damping of compact star oscillations 
previously took into account only the first, sub-thermal 
term in the parameterization eq. ( 45 1 . The simple ana- 



lytic form allows one to conveniently extend these stud- 
ies in order to include large amplitude effects encoded 
in the second term. The small deviations of the sim- 



plified parameterization eq. ( 45 1 from the exact value of 



the bulk viscosity are negligible compared to the consid- 
erable uncertainties inherent in such a damping analysis. 
Evaluation of this expression requires knowledge of the 
susceptibilities B and C that depend on the equation of 
state. 



C. Models of quark matter 

We now apply the results derived above to some sim- 
ple models of quark matter. We start with the simplest 
model, free quarks in a "confining bag". We will call this 
a "quark gas" (QG). We consider a 3-flavor quark and 
electron gas, with massless electron, up and down quarks 
and strange quark of mass m s with pressure 
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(46) 



where the strange quark Fermi momentum is given by 
/i 2 — to 2 . Here B is the phenomenological bag con- 



PFs 



stant that is important for the equilibrium composition 
of a strange star, but does not affect transport properties 
like the bulk viscosity studied in this work. The equilib- 
rium state is determined from eq. (46) by taking into 



account charge neutrality and weak equilibrium with re- 
spect to both the explicitly considered non-leptonic chan- 
nel as well as the quark Urea channel. 

In quark matter there are multiple channels for beta 
equilibration: as well as the nonleptonic channel (32) 
there are Urea channels which convert d or s quarks in 
to u quarks and electrons, and emit neutrinos. However, 
at temperatures and oscillation frequencies of interest for 
compact star physics the Urea rates are much slower, and 
their contribution to the bulk viscosity is heavily sup- 
pressed. This means that the fractions x u and x e remain 
constant during the oscillation. The required susceptibil- 
ities then are given by 
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Taking into account charge neutrality, the above equation 
of state yields to leading order in m s j '/i 9 the susceptibil- 
ities given in table [ill ( for the case c = 0) . 

According to eq. ( |43[ ) the maximum viscosity of a quark 
gas is given by 



Cr. 



12tt 2 w 



(49) 



which depends on density only through possible density- 
dependence of the strange quark mass. 

In Fig. [4] we give a comparison between the parameter- 
ization eq. (J45J) and the full numeric solution. We show 
the amplitude dependence of the viscosity of strange 
quark matter for a range of temperatures These results 
are analogous to those given by Madsen in his initial anal- 
ysis of supra-thermal effects |T7| . The analytic solution 
features the qualitative form that has been observed for 
the general result in fig.[2]and shows a striking agreement 
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Figure 4: The viscosity of a strange quark gas as a function 
of the amplitude of the density oscillation An/n for different 
temperatures. The plots are given for an intermediate density 
n = 2no and a frequency oj = 8.4 kHz corresponding to the 
oscillation frequency ui — 4/3 Q, of the quadrupole r-mode of 
a millisecond pulsar. The viscosity increases with the given 
temperatures starting from 10 6 K (bottom) to 10 9 K (top) 
and then decreases again. The thick, dashed curves represent 
the analytic model parametrization eq. J45} and the thin, full 
curves beneath them give the full numeric result. Clearly the 
parametrization is very accurate in the relevant regime below 
the maximum. At high temperatures the viscosity does not 
reach the supra-thermal regime for any physical value of the 
amplitude, hence the horizontal lines for T Js 10 10 K. 



with the full solution in the physically relevant region of 
amplitudes below the maximum. Note that for tempera- 
tures around T max the parametrization eq. ( |45| overesti- 
mates the viscosity for amplitudes above (An/h) max , as 
can be seen for the T — 10 9 K curve in fig. [I] However, 
if such amplitudes are reached then suprathermal bulk 
viscosity is overwhelmed, and other physics will have to 
be invoked to stop the growth of the mode. 

We now examine the sensitivity of our results to un- 
certainties in the quark matter equation of state. We use 
an extension of the phenomenological parameterization 
proposed in [ID] that allows us to study the behavior of 
the equation of state around chemical equilibrium. Ex- 
panding the ideal gas pressure to quartic order in m s , the 
m s -independent quartic terms in the individual quark are 
modified 
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Table II: Strong interaction parameters describing the response of various models of dense matter. In the case of hadronic 
matter with baryon density n a quadratic ansatz in the proton fraction x parameterized by the symmetry energy S eq. (1571) is 



employed. The expressions for a free hadron gas are given to leading order in n/m%, and for quark matter with quark chemical 



potential [i q using eq. ([50]) to next to leading order in m a /fj, q . The parameter c takes into account interaction effects within 
the employed quark matter model and vanishes for an ideal quark gas. 



where c is a new parameter which incorporates some ef- 
fects of strong interactions between the quarks and m s 
can parametrize here, in addition to corrections arising 
from the strange quark mass, also other interaction ef- 
fects, like the pairing gap in color superconducting mat- 
ter @Q]. 

The bulk viscosity is sensitive (via the susceptibilities) 
to the parameters c and m s , but not to the bag con- 
stant. We show in fig.[5]the effect on the bulk viscosity of 
varying c and m s within their expected range of values, 
at twice nuclear saturation density and a temperature 
T = 10 8 K. We calculate the bulk viscosity for an angular 
frequency of the oscillation of u = 8.4 kHz (correspond- 
ing to the r-mode of a pulsar with a period of 1 ms) . We 
find that the uncertainty amounts to more than an order 
of magnitude. In contrast to the equilibrium composition 
of strange stars which proved to be strongly dependent 
on the parameter c [40j, in the present case the effective 
strange quark mass has a larger impact. 

Finally we show in fig. [6] the dependence of the viscos- 
ity of a quark gas on the density of the matter and the 
frequency of the oscillation. The density dependence is 
most pronounced in the sub-thermal regime and becomes 
basically irrelevant in the supra-thermal regime, in accor- 
dance with the density-independence of the maximum of 
the viscosity eq. (49 1. Further, we see that the viscos- 
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ity increases strongly with frequency, according to the 
1/oj-dependence of the maximum eq. (49 I which arises 
as a prefactor in eq. |45|. Therefore, the results for a 200 MeV (highest two curves, cyan online) 



Figure 5: The dependence of the viscosity on parameters 
of the equation of state of strange quark matter using the 
simple parameterization eq. (50). We show the amplitude 
dependence at T — 10 8 K for lo — 8.4 kHz and n — 2rio. 
Dashed curves are for c = 0, solid curves are for c = 0.3. 
We show m 3 — 100 MeV (lowest two curves, magenta online) , 
m s = 150 MeV (middle two curves, blue online) and m s = 



millisecond-pulsar given here in all other figures present 
a lower limit for the viscosity, whereas the damping of 
slower rotating stars is much faster. 



IV. HADRONIC MATTER 

A. General features 

The bulk viscosity has been calculated for various 
phases of nuclear matter (unpaired, superfiuid, kaon- 



condensed etc) with flavor equilibration via either direct 
or modified Urea processes |19fB4] . The leptonic contri- 
bution has recently been calculated |41| . and hyperonic 
matter has also been studied [25H25] . We concentrate on 
the simplest case of non-superfluid hadronic npe matter. 
We note however, that the generic properties of our re- 
sults also apply to more complicated forms of matter like 
hyperonic and/or superfiuid nuclear matter. In the case 
of hadronic matter we assume that weak equilibration 
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Figure 6: The dependence of the viscosity of quark matter on 
the density and oscillation frequency, using the phenomeno- 
logical equation of state eq. (50 1 with m s = 150 MeV and 
c = 0, at T = 10 8 K. The dashed (blue online) "baseline" 
curve is for n = 2no and a high angular frequency cj = 8.4 
kHz corresponding to a millisecond pulsar. The dot-dashed 
(orange online) curves show the variation from the baseline 
with density: a low value n — 0.5no and a high value of 
n — 5no- The dotted (purple online) curves give the varia- 
tion from the baseline with angular frequency: a lower value 
uj — 0.84 kHz and an intermediate value uj = 2.8 kHz. 



occurs via the Urea channel 



P- 



p + e + v e 



(51) 



There are two qualitatively different cases depending on 
whether the direct process is possible or only the mod- 
ified version where a bystander nucleon is necessary to 
satisfy energy-momentum conservation. The latter rep- 
resents a particular strong interaction vertex correction 
to the above process. However, from the point of view 
of the weak interaction these different processes belong 
to the same channel. Taking into account baryon num- 
ber and charge conservation Sn p — 5n e = 0, the driving 
baryon number density oscillation yields here the oscil- 
lating chemical potential difference 



Pi = Mn — Pp ~ Pe 



(52) 



where the notation reflects that the equilibrating quan- 
tity in this case is isospin. 

Taking into account the effect of supra-thermal oscil- 
lation amplitudes requires the non-linear /i/-corrections 
to the corresponding rates. These have been given for 
hadronic matter in |191 W2\ |4"3"] in the standard case that 



only modified Urea processes are allowed 
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and in the enhanced case when direct Urea processes 
dominate HUH 
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where T 8 is the temperature in units of 10 8 K. Here 
we use the expressions given in |42j . but we note that 
the hadronic rates depend on model assumptions for the 
behavior of the strong interaction at high density (see 
[431 FH] ) and thereby involve uncertainties. These expres- 
sions yield the parameter values given in table IT] There 
are major differences between these hadronic rates and 
the corresponding one for strange quark matter. In quark 
matter, non-leptonic processes are naturally allowed and 
only particles that have a Fermi surface (quarks in this 
case) are involved. In contrast in hadronic matter such 
processes are absent (unless hyperons are present) and 
equilibration must proceed via semi-leptonic processes 
involving particles with no Fermi surface (neutrinos in 
this case) giving a much stronger temperature depen- 
dence. As noted before the simple analytic approxima- 
tion suitable for strange quark matter is not applicable 
here. Nevertheless we will see that many qualitative as- 
pects of that solution obtain in the general case. We 
note that although the prefactors of the non-linear terms 
decrease strongly as the power of Hi/T rises, it is not 
sufficient to neglect them since they enter the non-linear 
differential equation (21) where they dominate at suffi- 



ciently large amplitudes. 

According to eq. (15), the susceptibilities for hadronic 
matter are 
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(55) 
(56) 



Computing these quantities requires the equation of state 
of dense neutron matter. We will perform calculations 
using two model equations of state of nuclear matter. 
The first one is the "hadron gas", consisting of an electri- 
cally neutral beta-equilibrated mixture of free neutrons, 
protons, and electrons. The second one is "APR hadron 
matter", using the well-known model by Akmal, Pand- 
haripande and Ravenhall |45| which relies on a potential 
model that reproduces scattering data at nuclear densi- 
ties. As a low density extension of the APR data we 
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use |4"6l, Wf\ . In order to make it easy to apply our gen- 
eral results to other equations of state, we implement 
the APR equation of state using the simple parameteri- 
zation employed in [48] to approximate the dependence 
of the energy per particle on the proton fraction x by a 
quadratic form 



Other features of the plot do depend on the equilibra- 



E(n, x) = E s (n) + S{n)(l~2x) 2 



(57) 



where E s and S are the corresponding energy for sym- 
metric matter and the symmetry energy. We perform 
a global quartic fit to the APR prediction for symmet- 
ric and pure neutron matter E n which then yields the 
symmetry energy as 



S(n) = E n (n) - E s (n) 



(58) 



and the complete pressure including the electron contri- 
bution reads 

. . , fdE s (n) dS(n) , „ , 9 \ at 



\ dn 
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12tt 2 
(59) 
In the absence of oscillations the /3-equilibrium condition 
fij = yields the electron chemical potential as 



u e =4(1- 2x) S(n) (60) 

and the requirement of charge neutrality n p — n e al- 
lows us to determine the proton fraction x(n) so that the 
pressure becomes a function of the baryon density alone. 
With these explicit expressions the susceptibilities in ta- 
ble [TTJcan be computed and the general results in section 
ITT] can be employed. In the following subsections we will 
discuss the numerical results for the bulk viscosity of nu- 
clear matter, comparing it with those for one particular 
model of quark matter, the one given by eq. (50 1 with 
to, = 150 MeV and c = 0.3. 



B. Sub-thermal case 

When /ia «Twe obtain from the analytic expression 
eq. (26 1 the results shown in fig. [7] where the bulk vis- 
cosity of strange quark matter discussed in the previous 
section is also included for comparison. Here and in the 
following plots we study matter at twice nuclear satura- 
tion density, n — 2n . and a compression cycle with a 
high angular frequency u — 8.4 kHz corresponding to an 
r-mode in a pulsar with a period of 1 ms. We see in fig. [7] 
that the maximum bulk viscosity of hadronic matter as 
a function of temperature (or equivalently as a function 
of angular frequency) is roughly an order of magnitude 
smaller than the maximum value for strange quark mat- 
ter. This is unrelated to the beta-equilibration rate: the 
maximum viscosity depends according to eq. ( 27 1 on the 



tion rate. As we expect from (26), quark matter achieves 



its maximum viscosity at the lowest temperature, and 
has less suppression at low temperatures. This is because 
the nonleptonic equilibration only involves two particles 
in the initial and final state, each of which has a large 
Fermi momentum ~ fi q and hence large phase space fac- 
tors. This leads to a low k — 2 and a large value of T 
(table[l|. Thus the suppression at low temperature is only 
T 2 , and, according to eq. ( 27 1 , T max is relatively low. The 



next fastest is the direct Urea process in nuclear mat- 
ter, which involves more particles (including neutrinos 
which have no Fermi surface and thus very little phase 
space) and therefore has a higher k and lower T, giving 
it stronger T 4 suppression at low temperatures, and a 
higher T max . The slowest is the modified Urea process in 
nuclear matter, which involves additional spectator nu- 
cleons, raising k to 6 and further lowering T, raising T max , 
and increasing the low-T suppression to T e . 

Note that the right-most solid and dashed curves in 
fig. [7] for hadronic matter with modified Urea equilibra- 
tion, correspond roughly to the leftmost of the three solid 
(red online) curves in fig. u[a) that run along the surface 
from front to back. 

We draw two important conclusions from fig. 17] First, 
we have retained the full resonant structure of the vis- 
cosity compared to previous analyses |19l H4"] where a 
low temperature approximation TBT K <C to was used. 
This allows us to see that the viscosity decreases again 
at large temperatures and the maximum (27 1 occurs at 



millisecond-scale frequencies at potentially physically rel- 
evant temperatures of the order 10 10 K for direct Urea 
and 10 11 K for modified Urea. This means that the res- 
onant structure may be important in some astrophysical 
applications and from eq. ( 27 ) it is clear that it becomes 



relevant susceptibilities of the matter in question. 



increasingly important at lower frequencies. Second, we 
see in fig. [7]that for nuclear matter there is a considerable 
difference between the solid curves which are based on 
an interacting equation of state |45H47| and the dashed 
curves which are for a free gas 3 of nucleons and elec- 
trons. These models have different susceptibilities B and 
C, and the main effect of this is a vertical shift of the 
whole curve. The shift in T max is smaller because of the 
square root in eq. ( |27| . Hadronic matter with interac- 
tions has been considered (with a more simplified equa- 
tion of state) in [20 1 44] but many analyses |lll I4"9l I5U] 
rely on the simple analytic result given by Sawyer |19| 
which is based on the free gas expression. We see that 
these differ by roughly a factor of three for the given den- 
sity of n — 2tiq, but according to fig. [3] this difference 



3 Note that strictly speaking there are no modified Urea processes 
in an ideal hadron gas. Yet, for comparison with previous studies 
we use here the interacting matter expression for the rate but the 
ideal gas expressions for the strong susceptibilities. 

4 Note that the numerical prefactor given in [19] is too large by 
two orders of magnitude, see also |51| . 
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Figure 7: The sub-thermal approximation to the viscosity at 
vanishing amplitude as a function of temperature for u — 8.4 
kHz and n — 2no. The right (red) curves represent the stan- 
dard case of hadronic matter with modified Urea processes, 
the middle (green) curves hadronic matter when direct Urea 
process are allowed, and the left (blue) curves strange quark 
matter with non-leptonic processes. The dashed curves are 
for the free hadron and free quark models; the solid curves 
are for APR hadron matter, and interacting quarks eq. (501 
with m a = 150 MeV and c = 0.3. With APR nuclear matter 
the bulk viscosity is ~ 3 times larger than for the free hadron 
gas used e.g. in [19, 50ll51| . 



can increase strongly both at lower and higher density. 



C. The supra-thermal regime 

Beyond the sub-thermal limit, a numeric evaluation of 
eqs. (21 1 and (23 1 is required. As discussed in sect, [jj] 



the temperature and amplitude dependence of the bulk 
viscosity for a given form of matter can be expressed in 
terms of the function T(d, /) which was plotted in a form 
that is independent of the equation of state of hadronic 
matter with modified Urea processes in fig. [2] 

Using this result, we show in fig. [5] plots of the ampli- 
tude dependence of the bulk viscosity at two tempera- 
tures (left panel: T = 10 6 K; right panel: T = 10 9 K) for 
the various forms of hadronic and quark matter consid- 
ered in this paper. Here solid lines again show the results 
for interacting matter whereas the dashed lines show the 
free hadron/quark gas results. 

Note that the right-most curves in [8] for hadronic 
matter with modified Urea equilibration, corresponds 
roughly to the foremost of the three dashed (blue on- 
line) curves in fig. u[a) that run along the surface from 



left to right. 

At the lower temperature the viscosity reaches the 
supra-thermal regime already for small amplitudes, 
whereas at the higher temperature the sub-thermal 
regime extends to large amplitudes, giving a flat 
amplitude-independent plateau at low amplitudes. The 
stronger non-linear feedback in the hadronic cases leads 
to a significantly steeper rise that correlates with the 
largest power in eq. (19 1. Interestingly, despite these dif- 



ferences the maximum value reached by varying the am- 
plitude is still roughly the same as the maximum value in 
the sub-thermal limit eq. (27 1, as has been analytically 
found in the case of strange quark matter. This is im- 
portant since it means that oscillations are approximately 
equally damped at all temperatures once the amplitude 
becomes sufficiently large. The maximum arises for am- 
plitudes of the order 0.01, 0.1 and 1 for strange quark 
matter and hadronic matter with direct and modified 
Urea, respectively. The supra-thermal enhancement of 
the bulk viscosity is so strong, particularly for hadronic 
matter, that it could well provide the main saturation 
mechanism for unstable r-modes, stopping their growth 
at amplitudes that are below the threshold for other com- 
peting saturation mechanisms (e.g. non-linear hydrody- 
namics) but large enough to allow spin-down of a neu- 
tron star via gravitational radiation on astrophysical time 
scales. 



V. CONCLUSIONS 

We have studied the bulk viscosity of dense matter in- 
cluding its non- linear behavior at large amplitudes. In 
particular we give a general solution for the bulk vis- 
cosity of degenerate matter that is valid for arbitrary 
equations of state and retains its full parameter depen- 
dence in sec. [TT] This allows one to include these supra- 
thermal effects in a systematic r-mode analysis. In the 
supra-thermal regime we give a general analytic result 
for strange quark matter with non-leptonic processes in 
sec. |III| We found that the free hadron gas model of nu- 
clear matter, used for example in [11| R"9"| 1ST)] to com- 
pute the susceptibilities that enter the viscosity, is not 
accurate even in the sub-thermal regime, and may sig- 
nificantly underestimate the viscosity. Moreover, we find 
that the standard low temperature (high frequency) ap- 
proximation is not applicable for temperatures around 
10 10 K and the full resonant form of the bulk viscosity is 
required. We confirm previous results for the amplitude- 
dependence of the bulk viscosity of strange quark matter 
[T7] and find that these supra-thermal effects are para- 
metrically even more important in nuclear matter, be- 
cause of higher-order non-linearities in the amplitude- 
dependence of the Urea rate. 

The most obvious application of our results is to the 
damping of unstable r-mode oscillations in neutron stars. 
As the amplitude of the mode enters the supra-thermal 
regime the viscosity will increase steeply above the sub- 
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Figure 8: Comparison of the bulk viscosity of the different forms of matter studied in this work as a function of the density 
oscillation amplitude An/n. The frequency is cj = 8.4 kHz, corresponding to an r-mode in a millisecond pulsar and n = 2n$. 
Left panel: low temperature T = 10 6 K; right panel: high temperature T — 10 9 K. The dashed curves are for the free hadron 
and free quark models; the solid curves are for APR hadron matter, and interacting quarks eq. (501 with m s = 150 MeV and 
c = 0.3. The bottom (red online) curves represent the standard case of hadronic matter with modified Urea processes, the 
middle (green online) curves are for hadronic matter when direct Urea process are allowed and the top (blue online) curves 
are for strange quark matter with non-leptonic processes. Our calculations are valid only for An/n <C 1, but we show their 
extrapolation to higher amplitudes in order to compare with the qualitative general structure of the solution in fig. [2] Note 
that this plot uses a high oscillation frequency and that the viscosity is even larger at smaller values. 



thermal result and can exceed it by many orders of mag- 
nitude, but eventually it reaches an upper bound that 
is completely independent of the particular weak damp- 
ing process and depends only on susceptibilities of the 
dense matter in question. The viscosity then decreases 
at even larger amplitudes. We conclude that if r-mode 
growth is not stopped by the supra-thermal bulk viscos- 
ity before this maximum is reached then other non-linear 
dynamic effects [L2HI5] will be required to stop it. We 
have already performed initial exploratory calculations 
of r-mode damping times, and these suggest that over 
a significant region of parameter space supra-thermal 
bulk viscosity is sufficient to saturate r-mode growth at 
a finite amplitude parameter a max < 1, as was previ- 
ously assumed [TT] . This topic will be discussed in more 
detail in a forthcoming publication. There are several 
other directions in which our research could be devel- 
oped. Other equations of state for quark matter could be 
studied, for example the perturbative equation of state 
|52| , and also other phases with different equilibration 
mechanisms. The same is true for the various equations 
of state and phases of hadronic matter. Our analysis 
was for the case of a single equilibration channel, and 
it would be interesting to extend it to multiple channels, 



which may be relevant to both hadronic and quark matter 
(see appendix A of [3D] and Ref. [35)) . In quark matter 
the non-Fermi liquid enhancement of the Urea rate |53| 
should further increase resonant effects. The application 
of our results to r-modes in neutron stars also raises inter- 
esting questions concerning the correct treatment of the 
crust [M] , and possible modification of the radial profile 
of the r-mode due to strong radial dependence of the bulk 
viscosity in layered stars such as hybrid stars. 



Finally we note that at low temperatures the 
suprathermal enhancement of bulk viscosity becomes 
large, and the amplitude threshold for entering the 
suprathermal regime becomes low. This may not be rel- 
evant to the damping of r-modes because they are also 
damped by shear viscosity which becomes large at low 
temperature. But for other modes of compact stars, such 
as monopole pulsations, shear viscosity will not play such 
a significant role, and suprathermal bulk viscosity might 
be the dominant source of damping if external perturba- 
tions make the amplitude large enough. This might be 
relevant to old, cold, accreting stars in binary systems. 
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